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Dogs have been essential to life in the Siberian Arctic for over 
9,500 y, and this tight link between people and dogs continues in 
Siberian communities. Although Arctic Siberian groups such as the 
Nenets received limited gene flow from neighboring groups, 
archaeological evidence suggests that metallurgy and new subsis- 
tence strategies emerged in Northwest Siberia around 2,000 y ago. It is 
unclear if the Siberian Arctic dog population was as continuous as the 
people of the region or if instead admixture occurred, possibly in 
relation to the influx of material culture from other parts of Eurasia. To 
address this question, we sequenced and analyzed the genomes of 20 
ancient and historical Siberian and Eurasian Steppe dogs. Our analyses 
indicate that while Siberian dogs were genetically homogenous 
between 9,500 to 7,000 y ago, later introduction of dogs from the 
Eurasian Steppe and Europe led to substantial admixture. This is clearly 
the case in the lamal-Nenets region (Northwestern Siberia) where dogs 
from the Iron Age period (~2,000 y ago) possess substantially less an- 
cestry related to European and Steppe dogs than dogs from the medi- 
eval period (~1,000 y ago). Combined with findings of nonlocal 
materials recovered from these archaeological sites, including glass 
beads and metal items, these results indicate that Northwest Siberian 
communities were connected to a larger trade network through which 
they acquired genetically distinctive dogs from other regions. These 
exchanges were part of a series of major societal changes, including 
the rise of large-scale reindeer pastoralism ~800 y ago. 
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arly archaeological and genomic evidence from Zhokhov 
Island in Arctic Siberia indicates that dogs belonging to a 
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distinct lineage were an essential component of life in the Arctic 
for over 9,500 y (1, 2). This tight link between people and dogs 
continues in Siberian communities such as the Koryaks, Itel’mens, 
Chukchi, and Nenets, where dogs continued to be used for 
hunting, herding, and sledding among other activities (3-5). Re- 
cent genomic data obtained from Samoyedic-speaking communi- 
ties such as Nenets and Selkups suggest that during the Holocene 
they received limited gene flow from neighboring groups, in- 
cluding Steppe pastoralists (6, 7). Given that humans and their 
dogs often migrate and interact in parallel (8), it is possible that 
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Significance 


The Siberian Arctic has witnessed numerous societal changes 
since the first known appearance of dogs in the region 
~10,000 years ago. These changes include the introduction of 
ironworking ~2,000 years ago and the emergence of reindeer 
pastoralism ~800 years ago. The analysis of 49 ancient dog 
genomes reveals that the ancestry of Arctic Siberia dogs shif- 
ted over the last 2,000 years due to an influx of dogs from the 
Eurasian Steppe and Europe. Combined with genomic data 
from humans and archaeological evidence, our results suggest 
that though the ancestry of human populations in Arctic 
Siberia did not change over this period, people there partici- 
pated in trade with distant communities that involved both 
dogs and material culture. 


Siberian dogs also experienced limited gene flow from other 
populations. 

In contrast to the human genomic evidence, linguistic and 
ethnographic data suggest more dynamic processes. Specifically, 
these data suggest that Samoyedic-speaking peoples of Northwest 
Siberia migrated from southern Siberia, or neighboring regions of 
southeast Europe, to the Arctic as recently as ~3,000 to 4,000 y ago 
(9-12). In addition, archaeological sites such as Ust’-Polui in 
Northwest Siberia show evidence of iron and bronze metallurgy and 
isolated finds such as glass beads that were likely introduced from 
the Steppe, Black Sea, or the Near East (13-15). The presence of 
this material culture suggests that these communities participated in 
broad-ranging trade networks (13-15). The proposed migrations 
and exchanges of materials and practices potentially also involved 
dogs, which could have led to admixture, improvement, and ulti- 
mately to the establishment of modern Siberian dog lineages such as 
the modern Samoyed breed. 

To assess whether the Northwest Siberian Arctic dogs pop- 
ulation was continuous, or was instead marked by admixture 
(possibly in relation to the influx of material culture from other 
parts of Eurasia), we sequenced 20 ancient and historical Siberian 
and Eurasian Steppe dogs ranging in age from 11,000 to 60 y ago 
and in genomic coverage between 0.1x and 11.1x (Dataset S1). 
We then analyzed these genomes alongside publicly available 


ancient (n = 29) and modern (n = 120) canids (Fig. 14 and 
Dataset S2). 


Results 


Evaluating Broad Ancestry Patterns in Siberian Dogs. We first 
assessed whether Siberian dogs possessed similar ancestry using 
both principal component analyses (PCA; SI Appendix, Fig. S1) 
and phylogenetic analyses (Fig. 1B and SI Appendix, Fig. S2). 
These analyses recapitulate findings of previous studies indicating 
that modern and ancient dogs can be broadly classified into three 
major groups: West Eurasian, East Asian, and Arctic/Americas (8, 
16). The West Eurasian lineage includes ancient Near Eastern and 
Levantine, modern African, ancient and modern European, and 
newly sequenced Bronze Age Steppe dogs (8, 16). The East Asian 
lineage consists of modern dogs from China, Vietnam, Island 
Southeast Asia, Dingoes, and New Guinea Singing Dogs. The 
Arctic/Americas lineage includes modern Arctic breeds such as 
Greenland Sled Dogs and Siberian Huskies, ancient American 
dogs, mid-Holocene dogs from Lake Baikal, a 9,500-y-old dog from 
Zhokhov Island, dogs from the Iamal-Nenets region, and historical 
dogs from across Siberia sequenced for this study (Fig. 1B). This 
indicates that the newly sequenced Siberian dogs in this study 
possess genetic ancestry that was continuous for at least 9,500 y. 

We then made use of TreeMix to assess admixture between 
these lineages. These analyses, based on 21 ancient and 14 
modern dog genomes, broadly recapitulate our neighbor-joining 
tree (Fig. 1B and SI Appendix, Figs. S2 and S3 A and B) and 
support the results of previous studies (8, 17) showing that the 
ancestry of European dogs resulted from admixture between 
Near Eastern and Arctic dog lineages as early as 10.9 thousand 
years ago (kya). Our expanded dataset indicated that the 
9,500-y-old Zhokhov Island sample is closest to a 6,000-y-old 
sample from Lake Baikal, rather than to ancient North Ameri- 
can dogs (18) (SI Appendix, Figs. S2 and S4A). This result is 
consistent with the split between ancient American and Siberian 
Arctic dogs and between divergent Siberian lineages being older 
than the calibrated radiocarbon date of the Zhokhov Island 
sample (~9.5 kya) and with an end-Pleistocene or very early 
Holocene introduction (and subsequent long-term isolation) of 
dogs into the Americas (18, 19). 

Analyses of a published and a newly sequenced genome from 
the Mesolithic site of Veretye I (~10.9 kya; Fig. 1) in Northeast 
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Fig. 1. (A) Map of ancient dogs included in the study with sample name and age (kya) with an Inset map of the lamal-Nenets region of Northwest Siberia. 


Data from samples represented by circles were generated in this study with a mean genome coverage between 0.1 and 19.9x; triangles represent publicly 
available ancient dogs. The colors of the data points represent the D-statistic value of the form D (black jackal, sample; Zhokhov, ASHQ01). Red-shifted colors 
show a closer affinity to ASHQO1 (ancient Near Eastern dog), and blue-shifted colors indicate a closer affinity to Zhokhov (ancient Arctic dog). (B) A TreeMix 
phylogeny with five migration edges that are indicated by gray dotted lines. Each population contains between one and three individuals (S/ Appendix, Table 
S1). The color of the branches correspond to average D-statistic value in A. Complete models with the outgroup and edge weights, as well as models with 
additional edges can be found in S/ Appendix, Fig. S3. 
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Europe recapitulate previous findings showing that these dogs 
possessed ancestry related to both Arctic (66% and 71%) and 
Western Eurasian (34% and 29%) lineages (8, 17) (Figs. 1 and 2, 
SI Appendix, Figs. S3 A and B and S4 A and D, and Dataset S3). 
This result suggests that akin to all other modern and ancient 
(post-Mesolithic) European dogs sequenced to date, Mesolithic 
dogs in Europe already possessed both Arctic and Western Eur- 
asian ancestry. The fact that the ancient Siberian dog from Zho- 
khov Island, dated to ~1,000 y after the Veretye dogs, possesses 
no Western Eurasian ancestry, however, indicates that Western 
dog ancestry had not reached the Siberian Arctic by 9.5 kya. 


Establishing the Ancestry of Bronze Age Steppe Dogs. The lack of 
Western Eurasian ancestry in the Zhokhov dog does not pre- 
clude the possibility of later influx of dogs into the Siberian 
Arctic from neighboring regions such as Europe, the Near East, 
or the Steppe. To test this hypothesis, we first established the 
ancestry of Steppe dogs from the Bronze Age [5 to 3 kya (20)], 
since their ancestry has not been characterized, and they represent 
a potential source of admixture in Siberian dogs. Analyses of one 
publicly available (Samaral, 0.7x) and one newly sequenced 
Bronze Age Steppe dog (Ishkinino, 1.4x) revealed that they form 
a genetically homogenous group that is closest to West Eurasian 
dogs (Fig. 14 and SI Appendix, Figs. S1 and S5). In fact, TreeMix 
analyses indicated that these dogs show strong affinities to ancient 
Near Eastern dogs (SI Appendix, Fig. S3 A and B), and admixture 
graph modeling suggested that Steppe dogs possessed additional 
Near Eastern ancestry greater than that found in ancient Euro- 
pean dogs (Fig. 2). In addition, D-statistics detected no evidence 
of admixture from other populations (including from Arctic dogs) 
in either of the two Steppe dogs since their most recent common 
ancestor (>3.8 kya) (SI Appendix, Fig. S5). 

Additional modeling using F4 ratios suggested that the ancestry 
of Steppe dogs can be modeled as ~40% Arctic (represented by 
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the Zhokhov Island dog) and ~60% Near Fast (represented by 
Iranian and Levantine ancient dogs), reference SI Appendix. 
Previous studies of ancient human genomes indicated that a 
population related to Chalcolithic Iranians contributed ~40% of 
the ancestry of Bronze Age Steppe populations likely as a result of 
expansion of farming from the Near Fast into the Steppe (21). 
Although we cannot rule out the hypothesis that Steppe dogs were 
introduced from Europe, our analyses suggest that, as in Europe 
(8, 22), the Neolithic expansion of farming from the Near East 
into the Steppe also involved the dispersal of dogs. The analyses of 
a more recent (0.7 kya) Steppe dog genome (Bolgar1) from the 
medieval city of Bolgar in Tatarstan (Russia) indicates that al- 
though European ancestry increased over time (Fig. 14 and SI 
Appendix, Figs. S2 and S4 A and D), the Near Fastern ancestry 
likely remained in dogs from the Steppe at least until the Middle 
Ages (Fig. 2 and SI Appendix, Fig. S6 B and C). 


Testing for Admixture in Siberian Dogs. In contrast to Bronze Age 
Steppe dogs, our analyses found that there are varying degrees of 
affinities (as revealed by D-statistics) among ancient and his- 
torical Siberian dogs to Near Eastern (Western Eurasian) and 
ancient Arctic ancestry (Fig. 1 and Datasets S1 and S4). In fact, 
our PCA indicated that all Siberian dogs sequenced in this study 
fall along a cline between West Eurasian and Arctic dogs (SI 
Appendix, Fig. S1), suggesting that Siberian dogs possessed dif- 
ferent levels of West Eurasian ancestry. This influx of non-Arctic 
ancestry into Siberian dogs could be the result of the introduc- 
tion of dogs with either Near Eastern/Steppe ancestry, European 
ancestry, or both. Testing these scenarios is difficult, since pre- 
vious studies have shown that “European” ancestry itself resulted 
from admixture between Arctic and Western Eurasian (Near 
Eastern) ancestries (8). 

In order to ascertain the different ancestry components of the 
dogs that admixed with the Siberian dogs dated between 10.9 and 
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0.06 kya, we defined six separate models (four of which are 
displayed in Fig. 2A) using the Admixture Graph R package (23) 
in which we used representatives of Near Eastern, European, 
and Arctic dogs, and iteratively tested how well each of the 
models fit the genomes of 22 Siberian dogs. We did not include 
East Asian dogs in this analysis for two reasons. First, archaeo- 
logical evidence suggests that Northwestern Siberian communi- 
ties primarily interacted with communities from the Steppe, 
Black Sea, and the Near East (13-15). Second, the results of our 
own analyses (SI Appendix, Figs. S3 A and B, S7, and S8 and 
Dataset S4) and previous studies (8) indicate that although 
Arctic dogs (represented by the ancient Zhokhov dog) may have 
contributed to East Asian dog ancestry, there is little evidence of 
gene flow in the opposite direction (SI Appendix). 

To represent each of the three ancestral lineages in the models, 
we chose two representatives from the Near East (TepeGhela, 5.8 
kya; ASHQ01, 2.3 kya), two representatives from Europe (New- 
grange, 4.8 kya; HXH, 7.0 kya), and one representative from the 
Arctic (Zhokhov, 9.5 kya). We then defined a set of models to 
characterize the ancestry of each Siberian dog as having derived 
from either of the following: 1) solely the Arctic dog lineage, 2) 
both Arctic and European dog lineages, 3) both Arctic and Near 
Eastern lineages, and 4) exclusively European and Near Eastern 
dogs only (Fig. 2). In all models, European ancestry was modeled 
as a mixture of Near East and Arctic (8). We then applied all four 
possible combinations of representative genomes and tested the 
goodness of fit with each Siberian/Steppe dog. 


Ancient Lake Baikal Dogs (~7 kya) Possess Limited Western Eurasian 
Ancestry. We first applied this method to assess whether the three 
Lake Baikal dogs (7.4, 7.0, and 6.9 kya) possessed any non-Arctic 
ancestry. Our analysis indicates that three of our four models fit 
the data equally well for the two oldest, lower coverage (Baikal2 
and Baikal3; ~0.3x) dogs (Fig. 2). In the case of the youngest and 
higher coverage (Baikall; 2.2x) Lake Baikal sample, however, we 
found that models involving admixture from Europe did fit slightly 
better (Fig. 2). Though it is tempting to interpret these results as 
evidence for the arrival of dogs possessing non-Arctic ancestry be- 
tween 7.4 and 6.9 ka, this finding could also result from the lower 
coverage of the two older genomes. To further query whether these 
dogs possessed a history of gene flow at all, we performed an ex- 
haustive model search using qpBrute (17), in which all three Lake 
Baikal dogs were analyzed as a single population (SI Appendix). The 
results suggest that only ~3% of more than 20,000 tested scenarios 
fit the data, and of these, the majority (~60%) involved no ad- 
mixture from a Western Eurasian source into any of the Lake 
Baikal dogs since their split with Zhokhov. 

Although more complex models with admixture from a Western 
source postdating the common ancestor of Lake Baikal and 
Zhokhov Island dogs did fit the data (also evident when consid- 
ering only the youngest dog [6.9 kya]; Fig. 2), F4 ratio tests sug- 
gested that this contribution was marginal (~9% of the ancestry 
in the Lake Baikal dogs). In addition, D-statistics analyses based 
on alignments to a red fox (Vulpes vulpes) reference genome 
(VulVul2.2) indicated a slight reference bias when compared to 
the canFam3.1 assembly (European Boxer dog) (SI Appendix, Fig. 
S10 and Dataset S4). However, regardless of the reference ge- 
nome, D-statistics for Baikall showed statistically significant sig- 
nals for gene flow from West Eurasian dogs, but this signal was not 
seen in the older Baikal dogs (Datasets S3 and S4). Combined, our 
results suggest that Lake Baikal dogs potentially received a limited 
influx of ancestry, around ~6.9 kya from Western Eurasian dogs. 
Greater depth of coverage from these samples is required to more 
confidently assess this scenario. If correct, however, this would 
suggest that the spread of dogs from West and Southern regions, 
into the Siberian Arctic, began between 9.5 kya (the age of the 
Zhokhov dog) and 6.9 kya. 
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Increasing Levels of Western Eurasian Ancestry in Siberian Dogs over 
the Last 2 kya. We then applied our admixture graph—based 
method to assess admixture from Western Eurasian sources into 
later Siberian dogs from contexts dating between 2.0 kya to the 
present day, following a period of an apparent absence of di- 
rectly dated dog remains in Siberia between ~6.0 kya to 4.3 kya 
(24). In contrast to the earlier Lake Baikal dogs, all of the best 
fitting models for these dogs involved admixture from a Western 
Eurasian source (Fig. 2), which suggests additional influx of 
European or Near Eastern-like ancestry into Siberia after ~6.9 
kya (the youngest Lake Baikal sample). Although in most cases 
models involving admixture from a European source fit best, 
some Iron Age (~2.0 kya) and medieval (1.1 to 0.85 kya) dogs 
obtained from the Iamal-Nenets region in Northwest Siberia 
(Fig. LA) potentially also possessed additional ancestry related to 
that of Near Eastern dogs (Fig. 2 and SI Appendix, Fig. S8). The 
signal for an additional contribution from the Near East to the 
medieval Iamal dogs was slightly weaker when using alignments 
to the red fox genome instead of the dog reference but is still 
present (compare SI Appendix, Fig. S6 B and C). 

We further explored this scenario through an exhaustive 
model search using qpBrute (17) in which the Iron Age and 
medieval Iamal dogs were treated as two populations (SI Ap- 
pendix, Fig. S11). Although the number and the topology of the 
fitting models varied slightly across the different runs (SI Ap- 
pendix, Table S3 and Fig. S11), all models showed different in- 
dependent admixture from European and/or Near Eastern dogs 
into medieval and Iron Age Iamal dogs. In line with nonlocal 
materials including glass beads and various metal items (13, 15) 
recovered from the same archaeological sites, our results indi- 
cate that the communities living in Iamal during the Iron Age 
and medieval periods were connected to a larger trade network 
through which they acquired dogs with non-Arctic ancestries (SI 
Appendix, Fig. S8). 


The Emergence of the Modern Siberian Dog Lineages. We then 
assessed whether the genome of two dogs sampled ~100 y ago 
(1927 CE) from a Nenets community on the Iamal peninsula was 
related to Jamal dogs from medieval and Iron Age periods. 
Identity-by-descent pairwise distances indicated that the closest 
ancient Arctic dog relative of these two historical Nenets dogs 
was a medieval Iamal dog from Ust’-Voikar (SI Appendix, Fig. 
S4 B and D) suggesting a degree of population continuity from 
the medieval period through to the early 20th century. Inter- 
estingly, this analysis indicated that these two 100-y-old Nenets 
dogs were also closely related to Samoyed dogs, a modern spitz 
breed introduced to the United Kingdom in the nineteenth cen- 
tury from Siberia and which became popular among polar ex- 
plorers such as Scott and Shackleton (25). This result indicates 
that although there have been multiple admixture episodes of 
Northwest Siberian dogs, their Arctic ancestry component survives 
in the modern Samoyed breed. Similarly, we found that modern 
Siberian Huskies share an affinity with historical East Siberian 
dogs and ancient Lake Baikal dogs (25) (Fig. 1 and SI Appendix, 
Figs. S4 C and F and S7). Together, these results indicate that 
several popular modern Arctic breeds maintained significant lev- 
els of ancestry from a lineage established prior to 9.5 kya in 
Siberia. Importantly, non-Siberian ancestry in Siberian dog breeds 
is not a modern phenomenon but has instead been an ongoing 
process for at least 2,000 y. 


Conclusions 


Our results indicate that Arctic dogs likely evolved in near iso- 
lation from other dog populations until at least the mid- 
Holocene (~7 kya). These ancient Arctic dogs likely inhabited 
a large region of Siberia from the New Siberian Islands to Lake 
Baikal. Over the last 7,000 y, however, the evolutionary history of 
Siberian dogs has been influenced by multiple introductions of 
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dogs from the Eurasian Steppe and Western Eurasia and little to 
no gene flow from East Asian dogs. Some of these introductions 
coincide with periods of major transformations within Northwest 
Siberian societies, including the introduction of metallurgy to the 
Arctic (15), the advent of the use of reindeer for transportation 
~2 kya, and the rise of reindeer pastoralism ~800 y ago (26). 
Altogether, this suggests that these profound transformations in 
Northwest Siberia were linked with the importation of material 
culture (including dogs) from neighboring regions through the es- 
tablishment of large-scale trade networks. This influx of genetic 
variability associated with dog lineages adapted to farming (Europe) 
and pastoralism (Steppe) potentially resulted in behavioral and 
morphological changes in Arctic dogs that may have facilitated the 
transition from foraging to pastoralism in the Siberian Arctic. The 
generation, analysis, and interpretation of additional ancient dog 
genomes within their archaeological context will help to address 
questions related to the role that dogs have played in the long-term 
human occupation of the Arctic, including the advent of novel 
survival strategies more generally, and more specifically, reindeer 
domestication and the transition to large-scale reindeer pastoralism. 


Materials and Methods 


Samples used for this study were obtained from archaeological contexts (n = 
19) or ethnographic collections held at museums (n = 10). Detailed de- 
scriptions of the samples used for this study, including site information and 
context information, are presented in S/ Appendix. DNA was extracted from 
each sample in laboratories dedicated to work on ancient or degraded DNA 
materials at the University of Copenhagen, Swedish Museum of Natural 
History, and Trinity College Dublin. All PCR reactions were conducted in 
separate facilities to prepare shotgun sequencing libraries; screening to 
determine endogenous DNA content for each sample was performed with 
Illumina sequencing technology, accordingly 20 of the 29 samples were se- 
lected for deep sequencing on Illumina and BGISeq platforms. Each sample 
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was aligned to the CanFam3.1 reference dog genome (27) using the 
Burrows-Wheeler Alignment Backtrack algorithm (BWA aln) (28, 29), sub- 
sequently pseudohaploid calling was performed on the samples and a panel 
of publicly available canid samples with ANGSD (30) to be used for down- 
stream analyses. Complete details of methods used for sampling, laboratory 
work, and computational analyses in this study are provided in S/ Appendix. 
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browser/home) under project ID PRJEB42261 which is publically available for 
download (31). 
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